Demon-free quantum Brownian motors 



o 
o 

(N 



L. Machura, 1,2 M. Kostur, 1 P. Hanggi, 1 P. Talkner, 1 and J. Luczka 2 

institute of Physics, University of Augsburg, Universitatsstrasse 1, D-86135 Augsburg, Germany 
^Institute of Physics, University of Silesia, P-40-007 Katowice, Poland 
(Dated: February 2, 2008) 

A quantum Smoluchowski equation is put forward that consistently describes thermal quantum 
states. In particular, it notably does not induce a violation of the second law of thermodynamics. 
This so modified kinetic equation is applied to study analytically directed quantum transport at 
strong friction in arbitrarily shaped ratchet potentials that are driven by nonthermal two-state noise. 
Depending on the mutual interplay of quantum tunneling and quantum reflection these quantum 
corrections can induce both, either a sizable enhancement or a suppression of transport. Moreover, 
the threshold for current reversals becomes markedly shifted due to such quantum fluctuations. 
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Brownian motors are small physical machines that op- 
erate far from thermal equilibrium by extracting energy 
fluctuations to generate work against external loads 
They present the physical analogue of biomolecular mo- 
tors that direct intracellular transport and control mo- 
tion and sensation in cells 0] . In contrast to these molec- 
ular bio-motors, however, the molecular sized physical 
engines necessitate - depending on the nature of parti- 
cles to be transported and their operating temperature - 
a description that accounts as well for quantum features 
such as tunneling and quantum reflection. For this class 
of quantum Brownian motors, recent theoretical stud- 
ies |3J, Li have predicted that the transport becomes dis- 
tinctly modified as compared to its classical counterpart. 
In particular, innate quantum effects such as tunneling 
induced current reversals, power-law like quantum diffu- 
sion transport laws, and quantum Brownian heat engines 
have been observed with recent, guiding experiments that 
involve either arrays of asymmetric quantum dots 5] or 
cell-arrays composed of different Josephson junctions Q. 

Because of the mutual interplay of quantum mechan- 
ics, dissipation and non-equilibrium driving, the theoret- 
ical description of such nonequilibrium, dissipative quan- 
tum Brownian motor devices is notoriously difficult. The 
present state of the art of the theory is characterized 
by specific restrictions such as e.g. an adiabatic driv- 
ing regime, a tight binding description, a semiclassical 
analysis, or combinations thereof 0, 0- As such, the 
study of quantum Brownian motors is far from being 
complete and there exists an urgent need of further de- 
velopments. The analytic study of quantum Brownian 
transport for arbitrarily shaped spatially periodic ratchet 
potentials presents such a challenge. This goal is ad- 
dressed here within the strong friction regime, where the 
underlying quantum dynamics can be modeled by a re- 
cently put forward, ingenious quantum generalization of 
Smoluchowski dynamics 0, Q • 

Classically, a system coupled to a thermal bath at tem- 
perature T is described in terms of Langevin equations or 
corresponding Fokker-Planck equations . For a Brown- 



ian particle this yields the Kramers equation which in the 
strong friction limit reduces to the Smoluchowski equa- 
tion. In quantum statistical physics the description of 
Brownian motion dynamics is distinctly more intricate; 
it has been worked out, however, in some detail within 
limited generality using e.g. the assumption of a linear 
bath dynamics or a weak coupling limit. For the latter 
case, quantum master equations, e.g. of Lindblad form, 
have been derived [1(1 . 

Quantum Smoluchowski dynamics- Recent work 
within the strong friction limit shows that quantum 
Brownian motion can be described by a generalized 
Smoluchowski equation that accounts for leading quan- 
tum corrections pi la] . For a particle of mass M moving 
in the potential V(x), Anker hold et al. proposed a quan- 
tum Smoluchowski equation (QSE) for the diagonal part 
of the density operator p(t), i.e. the rate of change of the 
probability density P(x,t) = (x\p(t)\x) in position space 
x assumes the form 71: 



jM^P(x,t) 



^V^(x)P(x,t) + ^D eS (x)P(x,t), 

(1) 

where 7 denotes friction. The effective potential reads 

V cS (x) = V(x) + (l/2)AV"(aO, (2) 

wherein the prime denotes the derivative with respect to 
the coordinate x. The prominent parameter 

A = O/71-M7) ln(n/3 7 /27r), = l/k B T, (3) 

describes quantum fluctuations in position space and ks 
is the Boltzmann constant. The effective diffusion coef- 
ficient reads Q 



D eS (x) = D Ank (x) = pr\\ + XpV"(x)]. 



(4) 



Note that eq. Q is valid whenever k B T <C foy. 

This so derived quantum-Smoluchowski equation ex- 
hibits, however, a disturbing shortcoming: In clear con- 
tradiction to the validity of the second law of thermo- 
dynamics, eq. Q yields for an arbitrary, asymmetric 
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periodic ratchet potential V(x) of period L at zero exter- 
nal bias a non-zero, stationary average velocity (v) = JL 
(or equivalently, a non- vanishing probability current J). 
This is so, because the expression for the current reads 



L 



1 - exp [*(£)] 



lM J { f dxD^ k (x) exp[-*(i)] C L du exp[vf( w )] ' 

(5) 

where ^(x) — J Q duVl s {u) / D ^ n ^(u) upon inspection is 
non-periodic with ^(i) 7^ 0. Thus, a finite stationary 
drift emerges; i.e. a Maxwell demon seemingly is at work 
at stationary, thermal equilibrium. 

Demon- free quantum-Smoluchowski dynamics - Next, 
we put forward a clear-cut modification of the above 
quantum-Smoluchowski equation which does not cause 
such a fake perpetual motion phenomenon. First, we 
observe from theory that the leading strong friction 
quantum correction involves the second order derivative 
of the potential V(x), see (@J. Following prior works 
0, 01 we shall consistently neglect (in the high friction 
limit) higher order contributions in A, which in fact would 
involve also higher order derivatives of the potential. 
This new, modified quantum Smoluchowski equation (M- 
QSE) is derived from the following set of construction cri- 
teria: We seek a new diffusion coefficient that (i), in lead- 
ing order reproduces the previous result in Refs. 0>H|]> 
and (ii), does not exhibit a Maxwell demon behavior, i.e. 
the modified dynamics yields in thermal equilibrium a 
vanishing probability current, and additionally (iii), the 
dynamics reproduces the correct thermal quantum posi- 
tion probability for strong friction [l^. The construc- 
tion criteria (ii) of zero flux together with the correct 
leading order result for the thermal position probabil- 
ity in (iii) then fixes the form of the diffusion function 
F[V"{x)] = a~ l [l - bV'(x)}- 1 uniquely. The two con- 
stants a and b read explicitly a = and b = A/3. 

Upon an expansion of F[V"(x)] into a series in A the 
two diffusion functions do coincide in first order with re- 
spect to the quantity e(x) — \\0V"(x)\ < 1, as required 
by the condition in (i). Therefore, this improved mod- 
ified quantum-Smoluchowski equation Q is given by a 
modified diffusion, reading 



D eS (x) = D mod (x) =0- 1 [l- X0V"(x) 



(6) 



Note that from a mathematical viewpoint our thermal 
M-QSE-dynamics assumes the form of a Pade-like, 11011- 
perturbative result in place of 0. The thermal quantum- 
Smoluchowski stochastic dynamics in this strong fric- 
tion limit is thus equivalent to classical Brownian dy- 
namics within the effective potential J5J) and the new, 
state-dependent diffusion coefficient given in JUJ). The 
corresponding (M-QSE) Langevin equation reads in the 
Ito-representation [9j 



where the dot denotes the time derivative and is 
(classical) Gaussian white noise of vanishing mean and 
correlation < £(i)£(s) >= S(t — s). The above scheme is 
close in spirit with the approximation method of colored 
noise driven dynamics in terms of corresponding effective 
Markovian processes |l3|. 

Quantum Brownian motor transport - A finite trans- 
port emerges when the system operates far from thermal 
equilibrium 0. In the present context, we investigate 
overdamped, quantum Brownian motors 0, E} with the 
quantum fluctuations characterized by the parameter A 
in © • To this aim, we complement the thermal quantum 
dynamics in eq. (0 with a slowly waggling nonthermal, 
deterministic or random force T)(t), i.e. 



1 Mx = -V^(x) + ^2 1 MD mod {x) + V(t) ■ (8) 
In dimcnsionlcss form we then obtain 



v = -Ks(v) + V22?mod(y) i(s) + V( S ), 



(9) 



where the position of the Brownian motor is scaled as 
y = x/L, time is re-scaled as s — t/ro, with the charac- 
teristic time scale reading To = M7L 2 / AV (the barrier 
height AV is the difference between the maximal and 
minimal values of V(x)). During this time span, a classi- 
cal, overdamped particle moves a distance of length L un- 
der the influence of the constant force AV/L. The effec- 
tive potential is W cS (y) = W(y) + (l/2)\ W"(y), where 
the re-scaled potential W{y) = V(x)/AV = W(y + 1) 
possesses a unit period and a unit barrier height. The 
dimensionless parameter Ao = A/L 2 describes quantum 
fluctuations over the characteristic length L. For exam- 
ple, the value Ao = 0.01 means that, roughly speaking, 
the difference between quantum and classical fluctuations 
of the position of the Brownian particle is significant over 
distances of the order s/AqL = 0.1L. The re-scaled diffu- 
sion function T> mo d{y) reads, 



V mod {y) = ^ 1 [l-\ p W"{y)Y 



(10) 



The dimensionless, inverse temperature 0o — AV/ksT 
is a ratio of the activation energy in the non-scaled po- 
tential and the thermal energy. The re-scaled Gaussian 
white noise is £(s) = (LJ AV)£(t) and the re-scaled, non- 
thermal force reads fj(s) — (L/AV)n(t). 

As a specific realization, we next consider nonther- 
mal fluctuations modeled by Markovian, two-state noise, 
7y(s) = {—a, a}, that switches with a rate v between 
the levels a and —a. This problem can be solved ana- 
lytically in the adiabatic limit, i.e. if v — > 0. In this 
limit the stationary averaged dimensionless velocity reads 
(y) = J = (1/2) [J(a) + J(-a)], where 



J(a) 



1 - exp(-/3 a) 



1 Mx = -V^{x) + s /2 1 MD Taod {x) 



(7) 



1 v+i 
Jdy' D mld(y) e M-Po^(y,a,)} J dz exp\j3 ^(z,a)} 
y 

(11) 
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and 

*(y,o) = W(y) + (l/2)X W"(y) - (l/2)X f3 [W'(y)} 2 
- (l/A)X 2 p [W"(y)} 2 + a\ p Q W'(y) ~ ay. (12) 

Its classical behavior, i.e. A = 0, has been studied in 
Refs. E3. 
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FIG. 1: Stationary velocity {y) vs. the two-state noise ampli- 
tude a for both, a strongly damped quantum Brownian motor 
(solid) and its classical counterpart (dashed). The ratchet 
potential of unit barrier height, W (y) — — W"o[sin(27rj/) + 
0.25 sin(47ry)], with Wo = 0.454..., is depicted with the inset. 
The theory (dotted) in Ref. yields a nonphysical (although 
small) quantum-Maxwell demon behavior at small noise am- 
plitudes a; further away from equilibrium (a > 0.5) the QSE 
and the M-QSE predictions practically coincide within line- 
thickness. The chosen dimensionless inverse temperature is 
/3o = 5. 

The influence of the quantum corrections is presented 
with Figs. 1-3. The role of quantum noise enters via 
two functions: The effective potential W c s(y) and the 
effective diffusion function D mo d(2/)- The quantum cor- 
rection to the potential depends logarithmically weakly 
on temperature. The crucial correction stems from the 
diffusion which increases as the temperature decreases. 
The prominent quantum effects appear for lower temper- 
atures. In Figs. 1-3, we take for the re-scaled quantum 
fluctuations a parameter value of Ao = 10 -4 ln(10 3 /3o)- 
This choice assures that the quantum Smoluchowski 
regime is fully valid down to low temperatures of the 
order /3o ~ 10. In Fig. 1 we depict the current vs. the 
dichotomic noise level a. We deduce that the quantum 
corrections reduce the absolute value of the current (the 
maximal absolute quantum correction is \XqPqW" (y)\ = 
0.202). Note that the current value approaches zero for a 
vanishing noise amplitude a — * (solid line) and, as well, 
for very large amplitude a — * oo. The modification of 
the diffusion coefficient turns out to be essential for small 
amplitudes a of the non-equilibrium two-state noise; this 
regime describes the near-equilibrium behavior with the 
directed current approaching zero. In clear contrast, the 



use of the conventional quantum Smoluchowski equation 
(QSE) in (dotted) yields a nonphysical, (although 
small) positive current value. It should be pointed out, 
however, that far away from equilibrium (for a > 0.5) 
the two forms of quantum-Smoluchowski dynamics yield 
practically identical results. 

The analytic expression for the current allows one to 
study arbitrarily shaped ratchet profiles. As an example 
we consider the more complex shaped asymmetric peri- 
odic potential 

W(y) = W {sin(2iry) + 0.4 sin[47r( 2 / - 0.45)] 

+0.3sin[67r(?/- 0.45)]}, (13) 

where Wo = 0.371 normalizes the barrier height to unity, 
see inset in Fig. 2. This ratchet potential exhibits an in- 



0.025 




0.1 1 10 100 

a 



FIG. 2: The dependence of the stationary current {y) vs. 
the two-state noise level a is depicted for the potential I13H 
(inset), for both, the modified quantum-Smoluchowski (M- 
QSE) theory (solid), the conventional quantum-Smoluchowski 
(QSE) theory (dotted) and the classical case (dashed), respec- 
tively, for an inverse dimensionless temperature /3o = 2. Note 
that the maximal absolute correction is actually rather small: 
\\ f3 W"(y)\ =0.06. 

triguing current reversal vs. the noise amplitude a. The 
maximal absolute quantum correction is |Ao/3oW"(?/)| = 
0.06. There occur two regimes: one regime of small noise 
levels a for which the amplitude of the quantum current 
is enhanced and one at larger noise amplitudes where 
the classical current exceeds its quantum counterpart. A 
most salient intermediate regime occurs for which the 
classical current is positive while the quantum current 
remains negative. The point of the physically relevant 
quantum current reversal is shifted towards larger noise 
levels. Use of the quantum-Smoluchowski diffusion the- 
ory of Refs. 0, @ yields a fake, positive-valued current 
at small dichotomic noise strength (dotted line), being 
accompanied by a non-physical (!) current reversal, see 
in Fig. 2. We find again that the expected convergence 
between the two theories occurs far away from thermal 
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FIG. 3: The directed quantum noise-induced transport (y) 
of the quantum Brownian motor (solid) vs. the dimensionless 
inverse temperature /3o for the ratchet potential see inset 
in Fig. 2, is compared with its classical limit (dashed) and the 
conventional quantum theory in J1J (dotted). The dichotomic 
noise level is set at a = 5. The variations of the quantum 
corrections are depicted with the dotted background. 



important advance over prior studies of quantum ratchets 
0,0] that often require the use of manageable, stylized 
potential forms. Our investigation additionally manifests 
a rich spectrum of quantum Brownian motor behaviors, 
exhibiting both, quantum induced enhancement and sup- 
pression of transport, as well as shifted current reversals. 

These novel features can advantageously be put 
to work for quantum ratchets on the micro- and 
nanoscale 0. Moreover, the structure of our quantum- 
Smoluchowski dynamics can be generalized to higher di- 
mensional overdamped situations as e.g. for quantum 
noise-induced directed transport on surfaces. In partic- 
ular, our method and these quantum ratchet signatures 
can be utilized to optimize transport properties in su- 
perconductors by controlling the motion of vortices and 
magnetic flux quanta |15L llfij | . 
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1517/13-4 and the ESF (Programme Stochdyn) is grate- 
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In Fig. 3 we elucidate the role of temperature. Equa- 
tion (|10fl shows that the quantum corrections increase 
monotonically as temperature decreases. We must re- 
frain ourselves, however, from analyzing the limit of ex- 
treme low temperature. This is so, because the quantum 
corrections then grow too large, causing the diffusion to 
pass from positive to nonphysical, negative values upon 
exceeding the threshold value 1; clearly, the strong fric- 
tion quantum theory is valid only below this threshold. 
In fact, for correction values close to threshold the non- 
diagonal, density matrix elements assume nonzero deco- 
herence values that can no longer be neglected with the 
quantum-Smoluchowski theory. Fig. 3 depicts these in- 
creasing quantum corrections with decreasing tempera- 
ture. At high temperatures the interplay between reflec- 
tion and tunneling transmission causes a larger (in ab- 
solute value) quantum current; upon crossing the point 
of classical current reversal this behavior is interchanged. 
At even lower temperatures quantum corrections cause a 
smaller current value. 

Conclusions.- By use of a distinct modification of 
quantum-Smoluchowski theory we have developed a 
strong friction quantum approximation that is in agree- 
ment with both thermal equilibrium statistics and - 
above all - with the second law of thermodynamics. This 
so obtained, modified quantum theory can be applied for 
far from equilibrium transport where it facilitates closed 
form expressions (in terms of quadratures) for directed, 
quantum Brownian motor transport. Our tractable re- 
sults hold true away from the semiclassical limit and, ad- 
ditionally, can readily be applied to experimentally, arbi- 
trarily shaped ratchet profiles. Note that this presents an 
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